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1.  INTRODUCTION 


The  Vulnerability/Lethality  (V/L)  Process  Stnicture  (Walbeit  et  al.  1993),  as  viewed  in  Figure  1,  is 
the  framewoilc  in  which  the  Army  now  conducts  V/L  analyses.  The  basis  for  the  taxonomy  of  the  V/L 
spaces  comes  ht)m  the  recognition  that  V/L  analyses  pass  through  distinct  levels  of  information  in  a 
precise  order.  These  levels  are: 

Level  1  -  ThreatA'arget  Interaction;  Threat/Target  encounter  conditions. 

Level  2  -  Taiget  Component  Damage  States. 

Level  3  -  Target  Capability  States,  and 

Level  4  -  Target  Combat  Utility. 

The  mappings  by  which  one  passes  from  one  level  to  the  next  are  dependent  on  different  kinds  of 
information  at  each  level.  Fc*  instance,  proceeding  from  level  1  to  level  2,  or  the  O,^  mapping,  describes 
the  physics  of  the  threat/taiget  interacdon.  Proceeding  from  level  2  to  level  3  (the  €>2^  mapping)  requires 
the  use  of  engineering  measurables/perfonnance  models  that  can  be  described  by  fault  trees.  Basically, 
fault  trees  describe  how  target  components  interact,  in  an  engineering  sense,  to  achieve  a  capability.  The 
Degraded  States  Vulnerability  Methodology  (DSVM)  is  one  example  of  an  O23  mapping  (Kunkel  1994). 
In  the  case  of  aircraft  flight,  the  ctq)abilities  are  dependent  on  such  environmental  factors  as  air  density 
and  mission  configuration.  Therefore,  to  implement  the  O2J  in  general  mission  independent  terms,  the 
C!q)abilities  must  be  expressible  in  the  form  of  functions  relaUng  the  capabilities  parametrically  to  these 
environmental  factors.  The  methodology  HELicopter  Performance  Evaluadon  (HELPE)  presented  in  this 
report,  provides  such  functional  reladonships  and  form  the  theoretical  foundation  for  aircraft  DSVM.  The 
logic  within  HELPE  and  DSVM  will  assist  in  the  development  and  evaluation  of  the  time-dependent  fault 
trees  that  will  be  incorporated  into  each  analysis.  Time-dependent  fault  trees  are  the  evaluation  of  damage 
mechanisms  such  as  fluid  leakers  (fuel,  lubrication,  and  hydraulic)  hrom  tanks  and  lines,  crack  propagation 
of  a  component  or  contamination  of  a  component  by  means  of  nuclear,  biological,  and  chemical  munitions. 
Eventually,  HELPE  and  DSVM  will  be  incorporated  into  the  new  stochastic  V/L  now  known  as  the 
Modular  Air-System  Vulnerability  Estimation  Network  (MAVEN). 

MAVEN  is  a  stochastic,  point  burst  model  (method)  for  rotary/fixed  wing  aircraft  and  missile  systems 
enable  of  both  vu  -rability/lethality  and  battle  damage  repair  (BDR)  analyses.  The  model  is  applicable 
during  all  [diases  of  the  acquisition  cycle  and  will  be  used  in  determining  predictions  for  live-fite  tests. 
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Schematic  reptesentations  of  the  vulnerability  wocess.  its  levels,  and  mappings. 
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Finally.  MAVEN  will  fill  a  void  in  the  current  vyL  analyses  by  providing  results  that  are  observable 
andAtr  tneasureaUe  (Roach  1994). 

Helicopters  are  widely  used  by  the  various  aimed  services  to  peifonn  many  of  die  unique  missions 
to  vdiich  they  are  adtqited.  Their  ability  to  take  off  and  land  vertically,  to  hover,  to  fly  dose  to  the  ground 
surface,  and  to  land  almost  anywiwre  make  them  ideal  to  perfoim  many  of  the  functions  diat  cannot  be 
perfbimed  by  other  aircrafL  In  dre  military  seivi(%s,  helicopters  are  used  for  close  support  of  ground 
oimbat  forces  and  to  quickly  transport  such  forces.  Helicopters  are  also  used  for  various  rescue  missions 
and  in  naval  operations.  Basically,  a  helicopter  is  an  aircraft  diat  uses  rotating  wings  to  generate  lift  to 
support  its  wei^  and  power  to  propel  it  throu^  the  air  and  to  control  its  flight 

Helicopters  are  designed  to  perform  specific  mission  objectives.  During  the  design  stage,  it  is 
desirable  to  be  able  to  evaluate  the  performance  capabilities  of  a  proposed  design.  For  threat  helicopters, 
one  also  wishes  to  assess  their  performance  ctqiabilities  on  the  basis  of  their  known  physical 
characteristics.  There  are  two  known  methods  for  evaluating  the  performance  of  a  helicopter.  The  first 
method  is  called  die  "force  balance  method"  and  die  second  method  is  called  the  "energy  balance  medrod." 
In  the  force  balance  method,  the  forces  on  the  different  blade  elements  of  the  rotors  are  integrated  to 
evaluate  the  net  rotor  force  and  the  net  torque.  The  helicopter  is  trimmed  to  force  and  moment 
equiliteium  to  maintain  a  qiedfied  fli^  condition.  The  determination  of  these  forces  and  moments 
requires  knowledge  of  die  rotor-induced  velocity,  as  well  as  the  rotor  blades  motion.  In  die  energy 
balance  method  of  analysis,  one  evaluates  the  power  requited  to  balance  the  individual  sources  of  energy 
losses.  Thus,  power  is  required  to  overcome  die  profile  drag  and  the  induced  drag  of  die  rotor.  Power 
is  also  required  to  overcome  die  parasitic  drag  and  to  provide  power  for  the  helicopter  to  dimb  vertically. 
The  energy  method  is  the  simpler  method  to  use.  Bodi  mediods  give  approximate  results.  One  of  die 
main  sources  of  errors  is  the  difficulty  of  accurately  accounting  for  the  many  sources  of  energy  losses  such 
as  at  the  transmission  rotor  interface,  the  rotor  hub,  die  tail  rotor,  the  landing  gear,  and  the  pods. 

Ifolicopler  performance  assessment  usually  addresses  the  following  questimis: 

•  For  a  given  helicopter  at  given  flight  conditions  and  specified  mission  wei^  how  fast  can  the 
helicopter  craise  in  level  flisbt  conditicms? 

•  Cm  tt  hover  in  die  out-of-ground  effea  mode  (OGE)  and/or  in  die  in-ground  effect  (IGE)  mode? 
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•  What  is  its  maximum  vertical  rate  of  climb  horn  a  hover  position? 

•  What  is  the  speed  for  maximum  endurance  in  the  air? 

•  At  what  speed  must  the  helicopter  fly  for  a  maximum  range? 

•  What  is  its  autoiotational  descent  ci^ability? 

The  computer  model  HELRE  addresses  all  the  questions  posed  above  except  for  die  autototadonal 
descent  capability. 

Implmnentiiig  this  model  within  DSVM  will  yield  capabilities  parametrically  as  a  function  of  energy 
loss.  exan^de,  die  parameters  include  such  items  as  net  rotor  force  and  net  torque.  The  questions 
above  will  satisfy  the  O23  m^iinng  as  viewed  in  the  V/L  process  structure. 

2.  THE  HELPE  COMPUTATIONAL  MODEL  SETUP 

HELPE  is  organized  in  a  main  program  (called  MAIN)  and  seven  subroutines  widi  descriptive  names. 
The  program  MAIN  plays  die  primary  role  in  the  computation.  It  reads  die  irqiut  and  writes  it  to  output, 
sets  up  the  omnputation,  and  calls  the  various  subroutines  needed  to  perform  different  aspects  of  die 
performanoe  evaluation.  The  {»indpal  part  of  the  computation  is  performed  in  diree  nested  DO  loops. 
The  first  lotqi  is  indexed  to  oondua  die  computation  for  die  first  NWT  cases  of  die  desired  qierational 
missitm  weirds.  A  one-dimensional  weight  array  (WT)  with  six  elements  is  defined  before  entering  the 
DO  loop  together  with  various  other  parameters  needed  to  define  the  (^lerational  crmditions  of  the  flight. 
A  one-dimensional  speed  array  VF  consisting  of  16  dements  r^Ncsenting  forward  levd  speeds  is  dien 
defined,  uliidi  covers  the  whole  raitge  of  possible  flight  speeds.  The  six  demmits  of  WT  are  defined  as 
the  overloaded  gross  weight,  die  gross  wdght,  the  average  of  the  gross  and  empty  wdghts,  and  die  oae- 
and  two-diirds  kiad  wdghts.  The  dements  of  VF  start  at  30  knots  and  increase  in  stqis  of  10  knots.  The 
subroutine  PAVAIL  is  then  cdled  to  compute  four  arrays  of  available  power  with  elements  corresponding 
to  the  forward  levd  speeds  array.  The  four  arrays  correspond  to  the  power  availatde  from  one  or  two 
entries  at  die  continuoos  power  (MCP)  or  at  die  intermediate  rated  power  (IRP)  settings.  A 

second  primary  DO  lot^  is  then  entered  in  which  the  power  rerpiired  for  each  of  NFV  (<=20)  levd  flight 
speeds  is  computed.  This  computation  is  basically  accomplisbed  by  calling  the  sutnoutine  CQINT.  This 
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subroutine  computes  for  a  given  coefficient  of  thrast  (CDi  and  advance  ratio  XMU  a  torque  coefficient 
CQO.  Subroutine  CX2INT  computes  the  torque  coefficiem  by  means  of  a  two-dimensional  inteipolation 
of  a  taUe  of  data  consisting  of  values  of  CQ  corresponding  to  values  of  XMU  and  CT.  This  table  of  data 
is  read  into  the  code  as  a  part  of  the  input  on  tiqw  unit  three  by  subroutine  RDATA.  Actually,  the  values 
of  Cr  in  the  table  are  multifdied  by  10*.  and  tiie  values  of  CQ  are  multiplied  by  10^.  This  table  of  data 
is  for  a  base  helicopter  system.  Two  corrections  are  then  applied  to  CQO  to  detennine  tiie  value  of  CQ 
for  the  candidate  helicopter.  The  first  correction  DCQ  accounts  for  tiie  additional  torque  required  because 
of  the  differences  in  the  frontal  drag  area  between  the  base  and  the  candidate  helicopter.  The  second 
correction  accounts  for  the  differences  in  the  efficiencies  of  tiie  rotor  system  in  generating  thnist  CQ  is 
then  multiplied  by  an  tqipropriate  scaling  power  factor  to  determine  the  power  required  for  the  horizontal 
flight  for  the  corresponding  speed  and  mission  wei^  Finally,  a  third  DO  loop  is  entered  in  which  the 
additional  power  required  for  vertical  climb  is  computed.  The  evaluation  of  tiiis  additional  power  is  made 
for  NCR  (<=10)  different  rates  of  climb  which  are  part  of  the  input  The  first  rate  of  dimb  should  always 
be  zero  and  should  correspond  to  level  flight  conditions.  Fbr  the  level  flight  case,  subroutine  CTINT  is 
called  to  compute  for  the  available  powers,  the  peifonnance  ciq>abilities  of  the  helicopter  and  its 
maneuverability.  CTINT  computes  the  maximam  thrast  that  can  be  generated,  the  load  factor,  tite  bank 
anj^e,  the  turn  radius,  the  rate  of  turn,  and  the  lateral  accderation  for  the  maximum  continuous  powo* 
available  and  the  maximum  available  power. 

Program  MAIN  also  evaluates  tiie  possible  craise  speeds  of  a  helicopter  cone^xmding  to  tiie  available 
power  settings.  Program  MAIN  also  computes  for  the  candidate  helicopter  the  power  required  to  hover 
and  to  do  a  vertical  dimb  from  theoretical  considerations.  This  tiworetically  computed  power  takes  into 
consideration  tiie  power  required  by  the  tail  rotor  to  trim  the  helicopter.  The  actual  calculations  ate 
described  later.  MAIN  finally  computes  the  hover,  maneuverability,  and  vertical  dimb  capability  in  tiie 
OCX  mode  with  two  engines  and  with  one  engine  power.  This  hover  and  vertical  dimb  performance  is 
then  repeated  in  tiie  ICE  mode.  This  cycle  of  calculations  is  repeated  for  eadi  of  tiie  NWT  mission 
weights  desired. 

3.  SUBROUTINE  FUNCTIONS 

This  section  describes  the  functions  performed  by  each  subroutine. 
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3.1  MAIN.  This  subfoutiiie  reads  the  input  defining  the  physical  characteristics  of  the  candidate 
helicopter  and  die  desired  operational  conditions  and  writes  it  to  output  It  calls  subroutine  RDATA  to 
read  and  write  tables  of  input  data.  It  sets  up  the  computational  flow  and  calls  subroutines  CQINT. 
PAVAOU  CnNT,  VINDRA.  and  GE  to  compute  some  aspects  of  the  performance  evaluations.  This 
module  also  evaluates  the  cruise  speeds  and  the  hover  and  vertical  dimb  ciqiabilities. 

3.2  RDATA.  This  subroutine  reads  three  tables  of  data  from  input  and  writes  them  to  output  The 
first  two  tables  consist  of  values  of  torque  coefficients  CQ  x  10^  as  functions  of  thrust  coefficients 
CT  X  10*  and  speed  ratio  XMU.  The  two  tables  corre^nd  to  two  different  RPM  settings  of  the  rotor 
for  the  OGE  mode.  The  third  table  is  for  hover  in  the  IGE  mode. 

3.3  PAVAIL.  This  subroutine  computes  the  power  available  as  a  function  of  speed.  The  four  arrays 
of  power  available,  PACl,  PAMl,  PAC2,  and  PAM2.  are  defined  for  the  available  power  fiom  one  and 
two  engines. 

3.4  CQINT.  This  subroudne  computes  the  power  required  by  die  helicopter  to  fly  at  a  certain  speed 
in  level  fli^  for  the  given  missitm  weight  bong  considered. 

3.5  CITNT.  This  subroutine  computes  the  level  flight  maneuver  ctqiability  of  die  helioopter  for  die 
availaUe  continuous  and  maximum  power. 

3.6  VINDRA.  This  subroutine  calculates  the  velocity  at  which  die  helicopter  win  have  to  fly  to 
remain  aiibome  the  maximum  amount  of  time.  This  is  the  velocity  that  corre^ionds  to  die  minimum  fuel 
flow  consumption.  This  subroutine  also  calculates  the  velociQ^  diat  wiU  give  maximum  flight  range. 

3.7  This  subroutine  calculates  for  die  hover  mode  in  IGE.  die  required  CQ  for  die  mis^on  CT. 
It  also  calculates  for  a  given  CQ.  dw  maximum  CT  that  can  be  generated. 

3.8  LACC.  This  subroutine  calculates  for  the  hover  mode  in  OGE,  the  maximum  dirust  that  can  be 
generated  for  the  power  available. 

3.9  GENLSO.  Tins  is  a  general  least  squares  fit  subroutine  of  data.  The  subroutine  is  taken  fiom 
die  system’s  math  library. 
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3.10  MATINV.  This  is  another  subioutiiie  from  the  system’s  math  library.  MATINV  is  used  by 
GENLSQ  for  inveiting  matrices. 

3.11  FNEOS.  This  is  another  subroutine  from  the  system’s  math  library.  FNEQS  is  used  by 
GENLSQ. 

4.  'THE  MATHEMA’nCAL  MODELS 

4.1  Mathematical  Models  and  Techniques.  Hus  section  discusses  some  of  the  madiematical  modeLc 
and  techniques  duu  were  used  in  this  code.  The  first  model  discussed  is  the  computation  of  die  availab 
engine  power.  Data  (see  Table  1)  of  the  power  output  of  die  T700-GE-7(X)  engine  at  various  altitudes  h. 
International  Standard  Atmosidiere  (ISA)  conditions,  and  zero  airspeed  indicate  diat  the  available  power 
ouqxit  is  a  linear  fimction  of  die  altitude  ’11"  This  data  shown  in  Table  1,  and  plotted  in  Hgure  2,  is 
derived  from  General  Electric  (1976). 


Table  1.  'r700-GE-700  Available  Engine  Power  as  a  Function  of  Altitude  ’’h” 


1  Altitude  h  (ft)  1 

I  ao 

5,000 

10,000 

15,000 

20,000 

IRP(HP) 

1 

1,360 

1,164 

960 

770 

MCP(HP) 

1  ^>254 

1,145 

1,000 

850 

720 

Sooioe:  Oaml  Electric  (1976X 


A  linear  regression  fit  of  this  data  ^ves  the  following  equations: 

IRP(h)  -  1546.8  -  0.0388  h  (1) 

MCP(h)  -  1266  -  0.0273  h  (2) 

Equations  1  and  2  are  only  qiplicable  to  the  T7(X)-GE'700  engine.  In  program  MAIN,  die  variables  HIRP 
and  HMCP  iqnesent  die  intermediate  rated  power  and  the  maximum  continuous  power  at  altitude  h  and 
are  eiqiressed  as  linear  functions  of  h.  The  sea  levd  values  of  these  variables.  SIRP  and  SMCP,  as  wen 
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T700-  CE- 700  ENGINE 
ISA  CONDITIONS 
0.0  AIR  SPEED 


Souice:  Oeaeril  Elednc  (1976). 

Hgure2.  Available  engine  power  output  as  a  fimction  of  altitude  "h." 

as  the  slopes  of  these  lines,  DIRPDH  and  DMCFDH,  are  Uieiefoie  pan  of  the  required  input  for  an  engine. 
The  ctmiputed  available  powers  at  altitude  h  are  assumed  to  be  at  ISA  temperatures;  hence,  a  conection 
most  be  made  if  die  temperatures  were  different  frmn  standard  ISA  conditions.  Again,  the  available  powo- 
is  assumed  to  be  a  linear  function  of  the  temperature;  hence;  the  slopes  of  die  functions  with  respea  to 
the  temperature  DIRPDT  and  DMCPDT  need  to  be  given  on  iiqiut  for  eadi  engine.  The  resulting 
avaOabte  power  functions  at  temperature  T,  for  die  engines  TIRP  and  TMCP,  are  dien  used  by  subroutine 
PAVAIL  to  define  the  availatde  power  at  each  speed  for  the  one  and  two  engine  operations  at  die  two 
power  settings  of  IRP  and  MCP.  Fbr  speeds  exceeding  SO  knots,  the  available  power  increases  slighdy. 
PAVAIL  uses  die  equatitms: 

PACl(I)  -  TMCP, 

PAMl(I)  -  TIRP,  for  VF(I)  ^  50.0  , 

PACl(I)  -  TMCP  +  0.42  (VF(I)  -  50.0),  and 

PAMl(I)  -  TIRP  +  0.62  (VF(I)  -  50.0)  For  VF(I)  >  50.0  .  (3) 
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These  equations  (3)  are  based  on  fits  of  the  data  for  the  T700-GE-700  engine  and  are  only  valid  for  it 

The  second  model  discussed  is  the  computation  of  die  power  required  for  a  givoi  level  flight  at 
speed  V  and  weight  W.  The  speed  is  represented  by  the  dimensionless  advance  rado. 


-  V/Vt. 

(4) 

2  11  R  .  N, 

(5) 

where  Vt  is  the  rotor  tip  velocity,  R  its  radius,  and  N  is  the  number  of  rotor  revolutions  per  second.  The 
initial  estimate  of  scaled  power  required  CQO  is  computed  by  subroutine  CQINT.  This  computation  is 
achieved  by  interpolating  p  and  CT.  This  computed  value  CQO  represents  the  power  required  1^  the  base 
helicopter  at  die  specified  CT  and  p.  Two  corrections  have  to  be  implemented  to  determine  the  power 
required  for  the  candidate  helicopter.  The  first  correction  is  needed  to  accoum  for  the  differences  in  the 
frontal  drag  area  between  the  base  helicopter  and  die  candidate  helicopter.  The  power  required  for  a 
helicopter  of  wei|d>t  W  to  fly  at  level  forward  speed  V  can  be  written  as  in  McCormidc  (1967): 

P  -  PO  +  OJ  wVCpAV)  +  0.5.p.f  •  (6) 

f  -  Cd  •  Af.  (7) 

The  first  term  on  die  right-hand  side  of  Equation  6,  Po,  represents  the  rotor  blades*  profile  drag;  die 
second  term  represents  the  induced  drag;  and  the  third  term  represents  the  parasitic  drag.  The  profile  drag 
varies  very  slowly  widi  but  can  be  considered  to  remain  iqiptoximately  a  constant  The  induced  drag 
term  dominates  at  low  speeds  but  becomes  negli^ble  in  comparison  wifo  the  parasitic  drag  at  higb  qieeds, 
or  equivaloidy  large  values  of  p.  Thus,  for  large  values  of  p  G.e.,  p  >  0.30), 

P  -  Pp  -  0.5p.f.  V^.  (8) 

Pp  is  die  parasitic  power  lost  to  drag,  p  is  the  air  density,  C^  is  die  body  drag  coefficient  and  Af  is  die 
fiontal  body  drag  area.  LetCQandCPbethecoefBdentsoftorqueandpowerdienfiomdidrusualdefiititions: 
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CQ  -  CP  -  P/  (iiR^pV,^). 


(9) 


where  R  is  the  main  rotor  radius.  Using  Equation  8.  Equation  9  reduces  to 

CQ  -  CP  -  OJ  Cj  Af  p^/dtR^).  (10) 

Differentiating  Equatitm  10  with  respect  to  while  holding  other  parameters  constant, 

ACQ  -  0.5  C^  Af  Ap^/CnR^); 

hence, 

Cj  -  2  icR2  ACQ/(Af  .  Ap^).  (11) 

Let  CQi*,  CQ2  corre^nd  to  the  large  values  pj  and  P2,  reiqrectively.  TTien 

C,  -  (2  hr’ /A,)  (CQj-CQij/fpJ-pj).  (12) 

Let  P2  =  0.36  and  Pi  a  0.34;  then  we  can  evaluate  C^  from  Equation  12  for  the  base  helicopter.  To 
calculate  the  correction  in  ttie  value  of  CQ  due  to  the  difference  in  the  frontal  drag  area  between  the 
candidate  helicopter  and  base  helicopter,  we  differentiate  Equation  10  widi  respect  to  the  ratio  A/(2  icR^) 
^^e  holding  p  constant  and  assume  that  C^  remains  constant  Thus, 

(ACQ)d  -  (CQ)c  -  (CQ)b  -  Ci  [(Af/2  K  R\  "  (Af/2  tcR^),,  ]  .  p^.  (13) 

The  c  and  b  subscripts  refer  to  the  candidate  ami  base  helicopter,  re^ctively.  Thus,  (ACQ)^  is  the 
correction  that  needs  to  be  applied  in  this  case.  (CQ),,  is  tiie  torque  coefficient  for  the  base  helicopter 
which  is  denoted  in  tire  code  by  CQO:  therefore. 
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(CQ),  -  (CQ)b  +  (ACQ)i. 


(14) 


The  second  conecticn  that  needs  to  be  made  is  to  accoum  for  the  differences  in  the  efficiency  of  the 
rotor  systems  of  die  candidate  and  base  helicopters  in  producing  thrust  This  is  accounted  for  by  dividing 
die  computed  (CX})^  by  an  efficiency  factor  e  (EFF  in  HELFE)  that  represents  the  efficiency  of  the 
candidate  hdicqiter  system  in  producing  thrust  relative  to  die  base  helict^r  system. 

CQ  -  (CQ),/e  -  [  (CQ)b  -  (ACQ)^  ]  /e.  (15) 

The  powers  required  for  the  horizontal  flight  denoted  by  P|j  is 

Ph  -  CQ  •  P.c*  (16) 

P„  -  *R^  .  p  .  /550.  (17) 

P„  is  the  power  scale  factor  and  is  given  in  horsepower.  The  power  needed  for  die  helicopter  to  dimb 
vertically  at  ftAi  is  given  by 


Pc  -  Vc  (W  +  OJ  Cj  .  A, .  p  .  v/)  /550.  (18) 

where  Af  is  the  top  flat  plate  area,  and  the  total  poorer  required  for  a  forward  dimb  is  approximately 

Pt  -  Ph  +  Pc-  (19) 

4.2  Level  Flight  and  Maneuverability.  For  die  case  of  levd  flight,  subroutine  CIINT  is  called  to 
calculate  die  maximum  thrust  diat  can  be  generated  and  die  ability  of  die  helicopter  to  turn  and  maneuver. 
Thus,  if  T  is  the  thrust  generated,  W  the  operating  wei^  Ij  the  load  factor  in  a  turn,  and  ^  is  the  bank 
angle  as  seen  in  Figure  3,  then  according  to  Saunders  (1975), 

If  -  T/W.  (20) 
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cos  #  -  W/T  -  1/1,. 


(21) 


tan 


(22) 


The  centrifugal  force  is  givoi  by 


Fe  -  (W/g)  .  (v2/R,); 


(23) 


hence. 


T  sin  ♦  -  (W/g)  (V^/Rj). 


(24) 


where  R^  is  the  turn  radius.  Equation  24  reduces  to 


Rj  -  /  (11.26  .  tan  ♦)  feet. 


(25) 


where  is  the  level  speed  in  knots.  Hie  rate  of  turn  0  in  degrees  per  second  is  found  to  be 


6  >  1091  .  tan  ^  /V^  deg/sec. 


(26) 


and  foe  horizontal  component  of  acceleration 


A|,  -  g  .  tan  #  .  ft/sec^ 
where  g  is  foe  gravitational  accderation. 


(27) 
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The  maximum  cruising  speed  is  found  by  determining  the  point  of  intersection  of  the  curve  of  the 
power  le^iired  vs.  speed  and  the  curve  of  the  power  available  vs.  speed,  as  shown  in  Hguie  4.  HELPE 
determines  the  maximum  possible  level  flight  speeds  corresponding  to  available  MCP  and  IRP  powers. 
Hgure  5  shows  gnydiically  the  method  of  computing  the  q)eed  of  maximum  endurance  in  air  V|,  and  die 
speed  of  maximum  range  Vj  corresponds  to  die  minimum  P,  on  the  power  required  curve  and 
corresporuls  to  the  point  of  minimum  P/V  on  the  same  curve. 


Hgure4.  Determination  of  maximum  level  flight  speeds. 
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POWER 


HgureS.  Determination  of  maximum  endurance  and  maximum  ranee  speeds. 


5.  INPUT  PARAMETERS 


The  following  parameten  define  the  required  input  to  tun  a  case.  Seven  iiqtut  cards  need  to  be 
defined  as  follows  with  their  formats. 


Card  1:  Fbnnat  (AlO,  3  I  10) 


TITLE  =  Assigned  name  of  aircraft  to  be  evaluated. 
NWT  s  Number  of  desired  mission  wei^its. 

NFV  s  Number  of  forward  speeds. 

NCR  =  Number  of  dimb  rates. 

NWT  <=  6.  NFV  <=  20,  NCR  <s  10. 

Card  2:  Foimat  (8  F  10.4) 

WTO  »  Gross  weigfit  of  helicopier  Qtounds). 

WTE  s  Empty  weigjit  of  helicopter  (pounds). 

OVLF  s  Overload  weight  factor  to  be  used. 

FA  s  Frontal  flat  plate  area  of  aircraft  (fee^). 

TFPA  s  Top  flat  ^ate  area  of  aircraft  (fe^). 
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SCFAB  =  Scaled  frontal  area  (FA/2jdi^)  for  base  aircraft. 

DTR  =  Distance  of  tail  rotor  shaft  from  main  rotor  shaft  (feet). 

Card  3:  Format  (8  F  10.4) 

ROTR  =  Main  rotor  radius  (feet). 

TRR  =  Tail  rotor  radius  (feet). 

TRTV  =  Tail  rotor  tip  speed  (ft/sec). 

TRA  =  Tail  rotor  cant  angle  (degrees). 

SIGMA  =  Main  rotor  solidity  =  ratio  of  blade’s  area  to  rotor  disc  area. 

SIGMATR  =  Solidity  of  tail  rotor  ratio  of  blade's  area  to  rotor  disc  area. 
ZNOENG  =  Number  of  engines. 

RPM  =  Main  rotor  rotational  speed  (nev/min). 

Card  4:  Fonnat  (8  F  10.3) 

ROTNO  =  Number  of  main  rotors. 

ALTl  =  Desired  altitude  for  performance  evaluation  (feet). 

TEMPI  =  Desired  temperature  for  performance  evaluation  (deg  F). 

TEMPS  s:  Standard  atmosphere  (ISA)  temperature  (deg  F)  at  ALTl. 

BNO  =  Number  of  blades  for  main  rotor. 

BNOTR  =  Number  of  blades  for  tail  rotor. 

EFF  =  Efficiency  factor  of  candidate  rotor  system  relative  to  base  rotor  system. 
DCLS  =  Parameter  representing  improvements  in  technology  of  blade  design. 

Card  5:  Format  (8  F  10.3) 

RC(i)  =  i*  rate  of  vertical  climb  (feet/min),  i=l,2,..NCR. 
where  RC(1)  is  recommended  to  be  taken  as  0.0. 

Card  6:  Format  (8  F  10.4) 

SMCP  =  Maximum  continuous  power  (HP)  of  engine  at  sea  level. 

SIRP  =  Intermediate  rated  power  of  engine  (HP)  at  sea  level. 

DMCTPDH  =  The  partial  derivative  9  (MCP)/9h. 

DIRPDH  =  The  partial  derivative  9  (IRP)/9h. 

DMCPDT  =  The  partial  derivative  9  (MCP)/9T. 

DIRPDT  =  The  partial  deri\  ative  9  (IRP)^. 

CD  =  Coefficient  of  drag  based  on  top  flat  plate  area. 

T  here  represents  the  atmosphenc  temperature. 

Card  7:  Format  (8  F  10.4) 

DTL  =  Ui^r  transmission  limit  (HP). 

XLTL  =  Lower  transmission  limit  (HP). 

FLAG  =  A  flag,  FLAG  =  0.0  for  base  helicopter, 

FLAG  *  0.0  for  any  other  helicopter. 

DO  s  Constant  term  in  the  quadratic  expression  for  profile  drag  coefficient  of  blade 
D1  =  Second  term  coefficient  of  in  the  quadratic  for  profile  drag  for  the  blade. 
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D2  «  Ooeflicient  of  in  the  quadratic  for  profile  drag  coefficient  for  blade. 

A  «  Slope  of  the  curve  for  coefficient  of  lift  vs.  angle  of  attack  for  the  main  rotor  blades. 

The  profile  drag  coefficient  is  Cdo  »  -t-  D1  .  ~  -t-  D2  .  where  «  is  the  angle  of  attack  expressed  in 

radians,  as  given  in  Gessow  and  Myers  (1967). 

Subroutine  RDATA  is  called  to  read  the  CQ'^IO^  vs.  CI^IO^  and  p  data  from  tape  input  The  read 
data  are  then  written  to  ouqwt  Samples  of  such  data  for  the  UH-60A  and  the  AH64  are  given  in 
Appendix  B.  The  iiqjut  data  are  read  as  follows: 

Card  1:  Foimat  (2  1 5) 

NMUD  s  number  of  p  data 
NCTD  s  number  of  CT  data 

Card  2:  Fbimat  (15  F  5.2) 

DMU(j)  =  The  j*  p  data,  j=l,...  J'lMUD. 

Card  3:  Fbnnat  (15  F  5.2) 

CTDO)  =  Hie  i*  CT  data.  ial.....JJCrD 

Card  4:  Foimat  (8  F  10.3) 

CQDlOj)  «  The  (ij)*.  CQ  data  of  set  1.  i=l,...J^MUD;  j=l....a^CrD. 

Card  5:  Foimat  (8  F  10.3) 

CQI>2(iJ)  *  The  (ij)*.  CQ  data  of  set  2,  i=l,....NMUD;  j=l....J^CrD. 

Card  6:  Fonnat  (2  F  10.2,  F  10.3) 

CITNOc)  =  the  k*  CT  for  IGE  data. 

CQIN(k)  =  the  k’*'  CQ  for  IGE  data. 

DINMU(k)  s  the  k***  MU  (p)  for  IGE  data 
k  s  1.2, — ,8. 

6.  APPUCATION  AND  COMPARISON  OF  RESULTS 

This  section  sets  up  die  computation  of  an  example  and  compares  the  computed  results  widi  available 
test  data.  Such  data  are  available  for  die  UH-60A  Hdicopter  udiich  is  shown  in  Hgure  6.  The  primary 
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Figure  6.  Schematic  drawini 


q)entianal  mission  weight  for  the  UH-60A  is  16^60  lb.  The  perfonnanoe  evaluatirai  is  desired  at 
4.000-ft  altitude  dx>ve  sea  level  and  at  the  hot  temperature  of  95**  F.  The  empty  wei^  was  taken  to  be 
13,620  lb.  Test  data  for  the  UH-60A  during  these  crmditions  can  be  found  in  Nagata  (et  aL  1981).  The 
ii^  card  images  used  in  this  case  were  as  follow: 

Card  1:  TITLE  =  UH-60A.  NWT  =  4.  NFV  =  18.  NCR  =  6, 

Card  2:  WTO  =  16260.0,  WTE  =  13620.0.  OVLF  =  1.245,  FA  =  1455. 

TFPA  =  206.0,  SCFAB  =  0.0165,  DTR  =  32J67. 

Card  3:  ROTR  =  26.83.  TRR  =  5.5,  TRTV  =  685.28,  TRA  =  20.0, 

SIGMA  =  0.0826,  SIGMATR  =  0.1875,  ZNOENG  =  2.0,  RPM  =  257.9, 

Card  4:  ROTNO  =  1.0,  ALTl  =  4000.0,  TEMPI  =  95.0,  TEMPS  =  44.7, 

BNO  =  4.0,  BNOTR  =  4.0,  EFF  =  1.0,  DCLS  =  0.0, 

Card  5:  RC(1)  =  0.0,  RC(2)  =  450.0,  RC(3)  *  600.0,  RC(4)  =  800.0, 

Ra5)  =  1000.0,  RC(6)  =  2000.0, 

Card  6:  SMCP  =  1254.a  SIRP  =  1540.  DMCPDH  =  -0.0265,  DIRPDH  =  -0.0387, 

DMCPDT  =  -3.9563.  DIRPDT  =  -3.6819,  CD  =  0.234, 

Card  7:  UTL  =  2828.0,  XLTL  =  2800.00,  FLAG  =  0.0,  DO  =  0.0087, 

D1  »  -0.0216,  D2  =  0.400,  A  =  5.9244. 

Hgure  7  shows  HELPS  computation  of  die  powm*  required  at  the  primary  misami  wei^  of  16,260  Ib 
for  level  flic^  as  a  function  of  die  speed.  Plots  of  the  maximum  continuous  powm*  (MCP)  and  the 
intermediate  rated  power  (IRP)  n^ien  both  engines  are  operating  are  dso  shown.  The  {dots  show  diat  the 
minimum  cruise  qieed  with  MCP  during  these  condititnis  is  19  knots  and  die  maximum  cruise  speed  is 
148  knots.  HELPE  computadtm  shows  this  speed  to  be  148.18  knots  and  the  maximum  possible  speed 
with  IRP  power  to  be  165.11  knots.  Test  results  show  die  first  speed  to  be  147  knots.  HELPE 
computation  diows  diat  the  speed  for  maximum  range  is  119.5  knots,  udiile  the  test  results  give  it  as 
133  knots.  HELPE  amputations  and  Hgure  6  show  that  die  IRP  single  engine  speed  is  113J  knots. 
D^e  test  results  place  it  at  112  knots.  HELPE  also  gives  die  MCP  maximum  endurance  speed  to  be 
80.6  knots,  but  no  test  data  are  available  for  amqiaristm.  Test  results  show  diat  the  UH-60A  cannot  fly 
at  any  speed  cm  a  single  eitgine  MCP  power.  HELPE  computation  shows  this  to  be  die  case,  and  Figure  7 
illustiates  dds  graphically.  HELPE  computadon  shows  that  the  maximum  rate  of  vertical  dimb  with  IRP 
power  at  4,000-ft  altitude  and  95**  F  temperature  is  572  fl/hiin  idiile  test  results  give  it  as  590  fifoiia 
Figure  8  fluws  plots  of  the  maximum  speed  vs.  operating  weight  as  computed  by  HELPE  whmi  operating 
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Figure  8.  Computed  madmum  speed  vs.  opcraiing  wdidit  for  the  UH-6QA  helicopter. 

wiA  MCP,  and  IRP  power  at  die  4,000-ft  aMtode  and  95^  F  temperature  condititnis.  Table  2  shows  the 
HELPE  computed  power  required  for  hover  and  vertical  dimb  as  a  function  of  the  climb  rate  and  for  tbe 
various  operadcmal  missitm  wei^its  considered.  Appendix  A  contains  the  FORTRAN  IV  listing  of  the 
(teveloped  code.  Appendix  B  contains  a  listing  of  torque  vs.  thrust  coefficients  data. 


Table  2.  Power  Required  for  die  UH>60A  Helicopter  for  Hover  and  Vertical  Qimb  (OGE) 


Altioide  B  4/XX)  ft.  Temperature  =  95'*  F 

Available  Power  MCP  »  1,898  HP,  IRP  »  2,400  HP 

Climb  rate 

WTOb) 

(ftAnin) 

13,620 

14,940 

16,260 

20,244 

0 

1,667 

1,892 

2,118 

2,851 

450 

1,852 

2,096 

2,340 

3,127 

1,914 

2,164 

3,219 

1,997 

2,255 

3,342 

2,080 

2,346 

3,465 

2,495 

2,801 

mESMm 

4,081 
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The  HEU%<a»nputed  results  compare  very  favorably  with  die  test  measured  data.  Most  of  the 
con^mted  results  are  within  10%  of  fli^  test  measured  data.  HELPE  is  easy  to  mn.  and  the  computation 
is  very  ecmimnicaL  A  consideratde  amoum  of  information  is  obtained  in  output,  whidi  has  not  been 
described  or  compared  here  because  of  the  lack  of  test  values.  The  code  can  be  particularly  useful  in  the 
dedgn  stage  for  parametric  study  or  in  threat  systems  for  threat  evaluation.  Future  code  improvements 
and  extensions  might  follow  in  the  future,  but  this  report  documents  the  currnit  status  of  the  code. 
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USnNG  OF  THE  FORTRAN  77  PROGRAM 
OF  THE  COMPUTER  CODE  HELPE 
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PROGRAM  MAIN 

DIMENSION  RC(IO) ,WT(6) , PRCLD (10) , PREMR (10) ,PRTTH(10) 

COMMON  DMU(15) ,CTD (12) , CQDl (14,7) , CQD2 (14,7) ,NMUD,NCTD 
COMMON  /BCT/XMU,PSFAC, SIGMA, VTIP, CT, CLSTL, VFTS, RHOl 
COMMON/PAVA/PAC1E,PAM1E,VF(20) ,PAC1 (20) ,PAM1 (20) , PAC2 (20) , PAM2 (20) 
COMMON  /BCT1/CTCP,CTMP, IX,PRHF(20) ,PCLF (20) ,PMLF(20) , EFF 
COMMON/IGE/CTIN(10) ,CQIN(10) ,XINMU(10) 

READ  INPUT  FROM  TAPES  OR  (INPUT  DATA  FILE)  AND  WRITE  IT 
FLAG=0.0  FOR  BASE  A/C, FLAG  NOT  0.0  OTHERWISE  , SCFA=FA/ (2 . 0*AD) 

SCFAB  IS  THE  SCALED  FRONTAL  AREA  FOR  THE  BASE  A/C (UH60  OR  AH64) 
PAC1E,AND  PAMIE  ARE  MAX  CONT  AVAILABLE  POWER, AND  MAXPOWER  FOR  ONE 
ENGINE  AT  ALT1,AND  TEMPI  CONDITIONS, AND  0.0  AIRSPEED 
DTR  IS  DISTANCE  OF  TAIL  ROTOR  FROM  MAIN  ROTOR  SHAFT. 

SLMCP,SLIRP  ARE  CONT.  AND  MAX  POWER  AT  SEA  LEVEL  AND  STANDARD  COND . 
DMCPDH,DIRPDH,DMCPDT,DIRPDT  ARE  THE  PARTIAL  DERIVATIVES  OF  THE  POWER 
WITH  RESPECT  TO  ALTITUDE  AND  TEMPERATURE. 

CHARACTER  TITLE  *10 

READ (*, 130) TITLE, NWT,NFV, NCR 

READ (*,113) WTG , WTE , OVLF , FA, TFPA, SCFAB , DTR 

READ (*, 113) ROTR, TRR, TRTV, TRA, SIGMA, SIGMATR, ZNOENG, RPM 

READ (*,110) ROTNO, ALTl , TEMPI , TEMPS,  BNO, BNOTR,  EFF,  DCLS 

READ(*,110)RC(1)  ,RC(2)  ,RC  (3)  ,RC  (4)  ,RC  (5)  ,RC  (6) 

READ (*,113) SLMCP , SLIRP , DMCPDH, DIRPDH, DMCPDT, DIRPDT, CD 
READ ( * , 1 1 3 ) UTL , XLTL , FLAG , DO , D1 , D2 , A 
100  FORMAT (//) 

110  FORMAT (8F1 0.3) 

111  FORMAT (8F15. 4) 

112  FORMAT(2F10.2,2X,F10.2) 

113  FORMAT (8F10. 4) 

115  FORMAT (8F15. 4) 

120  FORMAT('CONTS  POWER  MAX  VEL=' , F8 . 2, lOX, ' MAX  POWER  MAX  VEL=',F8.2) 

130  FORMAT (AlO, 3110) 

140  FORMAT (' *******PERFORMANCE  CALCULATION  FOR  THE  " ,A8, '******' ) 

141  FORMAT ('AIRCRAFT=' , A8, 5X, 'NWT=' , I5,5X, 'NFV=' , I5,5X, 'NCR=' , IS) 

60  FORMAT ('ALT1=' , F8 . 0, 4X, ' TEMP1=' , F8 . 1, 4X, ' WEIGHT=' , F8 . 0) 

65  FORMAT (IH  , 6X, ' VHS' , 6X, ' CLM  RT',4X, 'REQ  HP',6X, 'AVAI  HP',7X, 'LO  FA 
*C',6X, 'BK  ANG',6X, 'TUR  R',6X, 'RT  TUR',4X, 'LAT  ACCEL') 

66  FORMAT (IH  , 5X, ' KNOTS' , 5X, ' FT/MIN' , 6X, ' HP' , lOX, ' HP' , lOX, ' RATIO'  ,  7X, 

*' DEGREES' ,5X, 'FEET' ,6X,' DEG/SEC' ,5X, 'FT/SEC2'  ) 

70  FORMAT (IH  , 1 OX, ' HOVER, VERTICAL  CLIMB  RATE,  AND  LAT  ACCEL  CAP  O.G.E 
*.  TWO  ENGINES') 

71  FORMAT (1 OX, 'HOVER, VERTICAL  CLIMB  RATE, AND  LAT  ACCEL  CAP  ONE  ENGINE 

*  O.G.E. ' ) 

72  FORMAT (1 OX, 'HOVER, VERT.  CLIMB  RATE, AND  LAT  ACCEL  CAP  I.G.E.,TWO  EN 
*GINES') 

73  FORMAT (1H0,10X,' HOVER, VERT. CLIMB  RATE, AND  LAT  ACCEL  CAP  I.G.E.,ONE 

*  ENGINE') 

75  FORMAT (IHO, 'AVAIL  CONT  POWER=' , E12 . 5, 3X, ' AVAIL  MAX  POWER=' , E12 . 5, 5 
*X, 'PREQ.  HOV.=' ,E12.5) 

80  FORMAT (IHO, 4X, 'CLIMB  RATE' , 5X, ' REQ.  PCLD. ' , 5X, ' REQ.  TOT.  POWER' , 3 
*X, 'REQ.  TH.  H.P. ' ,5X, 'LAT  CONT  ACCEL' , 5X, ' LAT  MAX  ACCEL' , 4X, ' CONT 
*VEL' ,4X, 'MAX  VEL') 
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81  FORMAT (IHO, 7X, 'FT/MIN' , 9X, ' H . P . ' , 13X, ' H . P . ' , 1 4X, ' H . P . " , 1 IX, 'FT/SEC 
*!• , 12X, 'FT/SEC2' , 9X, ' FT/MIN' , 5X, ' FT/MIN' ) 

84  FORMAT (IHO, 3E15.5, 3X,E15.5, 4X, E15 . 5, 5X, E15 . 5, 2X, F8 . 2, 5X,F8.2) 

82  FORMATdHO, 'CTCP=' ,E12.5,5X, 'CTMP=' ,E12.5,5X, 'CT=' ,E12.5, 5X, ' THRCP 

,E12.5,5X, 'THRMP=' ,E12.5) 

83  FORMAT (IHO, ' CTCI=' , E12 . 5, 5X,  '  CTMI  =  '  ,  E12 . 5,  5X,  '  CT='  , E12 . 5,  5X,  '  THRCP 
*1=' ,E12.5,5X, 'THRMPI=',E12.5) 

WRITE (*, 140) TITLE 

WRITE (*, 141) TITLE, NWT,NFV, NCR 

WRITE  (*,  111)  WTG,  WTE ,  OVLF ,  FA,  TFPA,  S CFAB ,  D TR 

WRITE (*, 115) ROTR, TRR, TRTV, TRA, SIGMA, SIGMATR, ZNOENG, RPM 

WRI TE ( * , 1 1 5 ) ROTNO , ALT 1 , TEMP 1 , TEMP  S , BNO , BNOTR, EFF , DCLS 

WRITE(*,115)RC(1)  ,RC(2)  ,RC(3)  ,RC(4)  ,RC(5)  ,RC(6) 

WRITE (*, 115) SLMCP, SLIRP, DMCPDH, DIRPDH, DMCPDT,  DIRPDT, CD 
WRITE (*, 115)UTL,XLTL,FLAG,D0,D1,D2, A 
C  SUBROUTINE  RDATA  IS  CALLED  TO  READ  CQ  VS  CT  DATA  FOR  CONSTANT  MU 

CALL  RDATA 

C  DEFINE  WEIGHTS  OF  INTEREST, PRESSURE, TEMP, DENSITY, THRUST  AND  POWER 
C  SCALE  FACTORS 

PI=3.1416 
AD=PI*ROTR*ROTR 
WT ( 1 ) =WTG*OVLF 
WT(2)=WTG 

WT(3)=0.5*  (WTG+WTE) 

WT(4)=WTE 

DW=WTG-WTE 

WT(5)=WTE  +  DW  /  3.0 
WT(6)=WTE  +  2.0  *  DW  /  3.0 
PRESS1«2117.- (0.1312*ALT1**0.9278) 

TRA1=459.7+TEMP1 
THETA=TRA1/518.7 
SCN=RPM/SQRT (THETA) 

RH01=PRESS1/ (32. 17*53. 34*TRA1) 

VTIP=6 . 283*ROTR*RPM/ 60 . 

TSFAC=3 .1416  *ROTR*ROTR*RH01 *VTIP  *VTIP  *ROTNO 
PSFAC=TSFAC*VTIP/550 . 

HIRP=SLIRP+DIRPDH*ALT1 
HMCP=SLMCP+DMCPDH*ALT1 
DT=TEMP1- TEMPS 
TIRP=HIRP+DIRPDT*DT 
TMCP=HMCP+DMCPDT*DT 
PAC1E=TMCP 
PAM1E=TIRP 

C  BEGINING  OF  THE  WEIGHTS  LOOP 
DO  600  IW=1,NWT 
WRITE (*, 100) 

WRITE (*, 60) ALTl, TEMPI, WT(IW) 

CT=WT (IW) /TSFAC 
CTOS=CT/SIGMA 

CALL  CQINT(0.34,CT,VTIP,CQ1) 

CALL  CQINT(0.36,CT,VTIP,CQ2) 

CDX=8948 . 35* (CQ2-CQ1) 

F=CDX*FA 
SCFA=FA/ (2.0* AD) 
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DSCFA=SCFA- SCFAB 
WRITE (*, 65) 

WRITE (*, 66) 

VF(1)=30.00 

C  BEGINING  OF  THE  HORIZONTAL  SPEEDS  LOOP 

DO  150  IV=1,NFV 
XIV=IV-1 

VF(IV)=VF(1)  +  XIV*10.0 
150  CONTINUE 

CALL  P AVAIL (NFV, ZNOENG, UTL, XLTL) 

DO  500  IV=1,NFV 
IX=IV 

VFTS=VF (IV) *1. 689 
XMU=VFTS/VTIP 
CLSTL=0.15-0. 07*XMU+DCLS 

C  CALCULATE  CQ  REQUIRED  FOR  HORIZONTAL  FLIGHT 

IF (FLAG  .EQ.  0.0) GOTO  200 
DCQ=CDX*DSCFA*XMU**3 
200  CALL  CQINT (XMU,CT,VTIP,CQ0) 

IF (FLAG  .EQ.  0.0)DCQ=0.0 
CQ= (CQO+DCQ) /EFF 
PRHF (IV) =CQ*PSFAC 

C  BEGINING  OF  THE  VERTICAL  CLIMBS  RATE  LOOP 
DO  400  IC=1,NCR 
RCS=RC(IC) /60.0 

PRCL=RCS* (WT(IW)+0.5*CD*TFPA*RHOl*RCS*RCS) /550.0 

PRTOT=PRHF (IV) +PRCL 

PRCLD(IC)=PRCL 

IF(RCS  .EQ.  0.0) GOTO  300 

WRITE (*, 112) VF (IV) ,RC(IC) ,PRTOT 

GOTO  400 

C  CALCULATE  THE  CTS  CORRESPONDING  TO  PAC  AND  PAM 
300  CALL  CTINT 
400  CONTINUE 
500  CONTINUE 

C  CALCULATION  OF  MAX  VELOCITIES  FOR  EACH  FLIGHT  WEIGHT 
DO  505  J=1,NFV 
IF(PCLF(J)  .GT.  1.0) GOTO  505 
IF(J  .LE.  2) GOTO  505 
GOTO  507 
505  CONTINUE 

507  SLC= (VF ( J) -VF (J-1) ) / (PRHF (J) -PRHF (J-1) ) 

VFCP=VF (J-1) +SLC* (PAC2 (J) -PRHF (J-1) ) 

DO  510  JJ=1,NFV 
IF(PMLF(JJ)  .GT.  1.0)  GOTO  510 
IF(JJ  .LE.  2) GOTO  510 
GOTO  512 
510  CONTINUE 

512  SLM= (VF (JJ) -VF (JJ-1) ) / (PRHF (JJ) -PRHF ( JJ-1) ) 

VFMP=VF (JJ-1) +SLM* (PAM2 ( JJ) -PRHF (JJ-1) ) 

WRI TE ( * , 1 2 0 ) VFCP , VFMP 
CALL  VINDRA (NFV, VF, PRHF) 

C  THEORETICAL  ESTIMATE  OF  POWER  REQUIRED  TO  HOVER  AND  TO  CLIMB 

B=1 . 0-SQRT (2 . 0*CT) /BNO 
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X=CT/ (B*B) 

CALCULATE  CQ  REQUIRED  BECAUSE  OF  PROFILE  DRAG 
CQPRO=SIGMA*DO/8. 0+0. 6667*Dl*X/A+4. 0*D2*X*X/ (SIGMA*A*A) 

C  CALCULATE  THE  CQ  REQUIRED  FOR  INDUCED  DRAG  AND  TO  CLIMB 

DO  513  J=1,NCR 
W=RC(J)  /60. 

VOVT=W/VTIP 

CQICL=0 . 5*CT*SQRT (VOVT*VOVT+2 . 0*X) +0 . 5*VOVT*CT 
CQDR=0 . 5*CD*TFPA*W**3/  (550 . 0*PSFAC) 

CQTMR=CQICL+CQPRO+CQDR 

THTR=CQTMR*TSFAC*ROTR/DTR 

TSCFTR=PI*RH01*'^RR*TRR*TRTV*TRTV 

CTTR=THTR/TSCFTR 

BTR=1 . 0-SQRT (2 . 0*CTTR) /BNOTR 

XTR=CTTR/ (BTR*BTR) 

CQPRTR=SIGMATR*D0/8 . 0+0 . 6667*Dl*XTR/A+4 . 0*D2*XTR*XTR/ (SIGMATR*A*A) 
CQITR=0.5*CTTR*SQRT (2. 0*XTR) 

CQTTR=CQPPTR+CQITR 
PRTTR=CQTTR*TSCFTR*TRTV/550 . 0 
PREMR ( J) *CQTMR*PSFAC 
513  PRTTH(J)=PREMR(J)+PRTTR 

C  END  OF  SECTION  OF  THEORETICAL  ESTIMATES  OF  POWER  REQUIRED 

C  BEGINNING  OF  HOVER, AND  VERTICAL  CLIMB  CALCULATIONS 

C  COMPUTATION  OF  POWER  REQUIRED  TO  HOVER  AND  TO  CLIMB  ON  BASIS  OF  DATA 

CALL  CQINT(0.0,CT,VTIP,CQHO) 

PRHO=CQHO*PSFAC/EFF 

C  NEXT  DO  LOOP  IS  FOR  FULL  AND  REDUCED  POWER  CALCULATION 
DO  550  K=l,2 
PAC*PAC1E*  ZNOENG 
PAM=PAM1 E * ZNOENG 
PAC=AMIN1 (PAC,UTL) 

PAM=AMIN1 (PAM, UTL) 

IF(K  .EQ.  2)PAC=AMIN1 (PAC1E,XLTL) 

IF(K  .EQ.  2)PAM=AMIN1 (PAM1E,XLTL) 

WRITE (*,100) 

IF(K  .EQ.  1)WRITE(*,70) 

IF(K  .EQ.  2)WRITE(*,71) 

WRITE (* , 75) PAC, PAM, PRHO 
CPCE5=PAC*EFF*1 . 0E+5/PSFAC 
CPME5=PAM*EFF*1 . 0E+5/PSFAC 
CALL  LACC(CPCE5,CPME5,VTIP,CTCP,CTMP) 

THRCP=CTCP  *  T  SFAC 
THRMP=CTMP * TSFAC 

WRITE (*, 82 ) CTCP , CTMP , CT , THRCP , THRMP 
WRITE (*,80) 

WRITE(*, 81) 

DO  522  J=1,NCR 
PRDAT=PRHO+PRCLD ( J) 

X1LF=CTCP/CT 

IF(X1LF  .LE.  1.0)ACCEL1=0.0 
IF(X1LF  .LE.  1.0) GOTO  515 
TANA1=SQRT (X1LF*X1LF-1 .0) 

DELPAC=PAC-PRDAT 
VCMAX=DELPAC*550 . 0/WT (IW) 
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VCL60=VCMAX*60.0 
ACCEL1=32 . 17*TANA1 

515  X2LF=CTMP/CT 

IF(X2LF  .LE.  1.0)ACCEL2=0.0 
IF(X2LF  .LE.  1.0) GOTO  520 
TANA2=SQRT (X2LF*X2LF-1 .0) 

DELPAM=PAM-PRDAT 
VMMAX=DELPAM*550 . 0/WT (IW) 

VML60=VMMAX*60.0 
ACCEL2=32 . 17*TANA2 

520  WRITE (*, 84)RC(J) ,PRCLD(J) , PRDAT, PRTTH ( J) , ACCELl, ACCEL2, VCL60, VML60 

C  NEXT  FOUR  STATEMENTS  ARE  TO  RETURN  TO  FULL  POWER  CALCULATIONS. 

522  CONTINUE 

550  CONTINUE 

C  IN  THIS  SECTION  WE  CALCULATE  THE  IN  GROUND  EFFECT  PERFORMANCE 

CALL  GE (CT, CQI ,0.0, CPCE5 , CPME5, CTCI , CTMI ) 

PRHOIG=CQI *PSFAC/EFF 
DO  590  1=1,2 
P AC=P AC IE*  ZNOENG 
PAM=PAM1 E  *  ZNOENG 
PAC=AMIN1 (PAC,UTL) 

PAM=AMIN1 (PAM,UTL) 

IF(I  .EQ.  2)PAC=AMIN1 (PAC1E,XLTL) 

IF (I  .EQ.  2)PAM=AMIN1 (PAM1E,XLTL) 

WRITE(*,100) 

CPCE5=PAC*EFF*1 . 0E+5/PSFAC 
CPME5=PAM*EFF*1 . 0E+5/PSFAC 
IF(I  .EQ.  1)WRITE(*,72) 

IF(I  .EQ.  2)WRITE(*,73) 

WRITE (*, 75) PAC, PAM, PRHOIG 

CALL  GE (CT, CQI, 1.0, CPCE5,CPME5, CTCI, CTMI) 

THRCP I =CTC I *  T SFAC 
THRMP I =CTMI * TSFAC 

WRITE (*,83)CTC I, CTMI, CT, THRCP I , THRMP I 
WRITE (*,80) 

WRITE (*,81) 

DO  580  J*1,NCR 
PRDIG=PRHOIG+PRCLD ( J) 

XLF1=CTCI/CT 

IF(XLF1  .LE.  1.0) ACCL1=0.0 
IF(XLF1  .LE.  1.0) GOTO  560 
TANA1=SQRT (XLF1*XLF1-1 .0) 

ACCL1=32.17*TANA1 
DPAC=PAC-PRDIG 
VCM=DPAC*550 . 0/WT (IW) 

VCM60=VCM*60.0 

560  XLF2=CTMI/CT 

IF(XLF2  .LE.  1.0)ACCL2=0.0 
IF(XLF2  .LE.  1.0) GOTO  580 
TANA2=SQRT (XLF2*XLF2-1 .0) 

ACCL2=32.17*TANA2 
DPAM=PAM-PRDIG 
VMM=DPAM*  550. 0/WT (IW) 

VMM60=VMM*60.0 
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580  WRITE (*, 84) RC(J) ,PRCLD(J) , PRDIG, PRTTH ( J) , ACCLl, ACCL2, VCM60, VMM60 
590  CONTINUE 
600  CONTINUE 
STOP 
END 

SUBROUTINE  RDATA 

COMMON  DMU(15) , CTD (12) ,CQD1 (14,7) ,CQD2 (14, 7) ,NMUD,NCTD 
THIS  SUBROUTINE  READS  THE  CQ  VS.  CT  DATA  FOR  A  GIVEN  VALUE 
OF  MU  FROM  TAPE  5  AND  WRITES  IT. 

COMMON/ IGE/  CTIN(IO) ,CQIN(10) ,XINMU(10) 

READ (*, 10)NMUD,NCTD 
READ (*,15)  (DMU(J) , J=1,NMUD) 

READ (*,15)  (CTD(I) , I=1,NCTD) 

WRITE (*,30) 

WRITE (*,11)  NMUD,NCTD 
WRITE (*,12)  (CTD(J) , J=1,NCTD) 

DO  101  I=1,NMUD 

READ (*,20)  (CQDl (I, J) , J=1,NCTD) 

101  WRITE(*,50)  DMU(I)  ,  (CQDl  (I,  J)  ,  J=1,NCTD) 

WRITE(*, 100) 

WRITE (*,40) 

WRITE (*,12)  (CTD(J) , J=1,NCTD) 

DO  102  I=1,NMUD 

READ (*,20)  (CQD2 (I, J) , J=1,NCTD) 

102  WRITE(*,50)  DMU(I) , (CQD2 (I, J) , J=1,NCTD) 

WRITE (*,60) 

DO  200  K*l,8 

READ(*,70)CTIN(K)  ,CQIN(K)  ,XINMU(K) 

200  WRITE (*, 75) CTIN(K) , CQIN (K) , XINMU (K) 

10  FORMAT (215) 

11  FORMAT (IHO, lOX, 'NMUD=' , 15, lOX, 'NCTD=' , 15) 

12  FORMAT (IHO, 8X,10F1 0.3) 

15  FORMAT (15F5. 2) 

20  FORMAT (8F1 0.3) 

30  FORMAT (IHO, 'OUT  OF  GROUND  EFFECT  TORQUE  COEFFICIENTS, CQ*1 . 0E5  VS  C 
*T*1.0E4' ) 

40  FORMAT (IHO, 'OUT  OF  GROUND  EFFECT  SECOND  SET  OF  TORQUE  COEFFS' ) 

50  FORMAT(1HO,F8.3, 10F10.3) 

60  FORMAT (IHO, ' IN  GROUND  EFFECTS  TORQUE  COEFF,  CT*1.0E4  VS  CQ*1.0E5') 
70  FORMAT (3F1 0.1) 

75  FORMAT(lH0,2F10.2,F10.3) 

100  FORMAT (//) 

RETURN 
END 

SUBROUTINE  CQINT (XMU, CT, VTIP, CQ) 

COMMON  DMU(15) , CTD (12) , CQDl (14,7) ,CQD2 (14,7) ,NMUD,NCTD 
C  THIS  SUBROUTINE  INTERPOLATES  FOR  THE  VALUE  OF  CQ  CORR.  TO  CT,MU 
CTE4=CT*1.0E4 
DO  30  I=1,NMUD 
IF (XMU  .GT.  DMU(I))GOTO  30 
GOTO  35 
30  CONTINUE 
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1=1-1 

35  DO  40  J=1,NCTD 

IF(CTE4  .GT.  CTD(J))GOTO  40 
GOTO  45 

40  CONTINUE 
J=J-1 

45  IF(J  .EQ.  1)J=J+1 

C  INTERPOL  FOR  CQ  ALONG  I=CONST  LINE (MU=COST) ;  J  AND  CT  VARY 

DCTJ1J=CTD ( J) -CTD (J-1) 

DCTJ1X=CTE4-CTD (J-1) 

SLOIl= (CQDl (I, J) -CQDl (I, J-1) ) /DCTJIJ 
CQI1=CQD1 (I, J-1)+SL0I1*DCTJ1X 
SLOI2= (CQD2 (I, J) -CQD2 (I, J-1) ) /DCTJIJ 
CQI2=CQD2 (I, J-1) +SL0I2*DCTJ1X 
IF (I  .EQ.  1) 1=1+1 

C  INTERPOL  FOR  CQ  ALONG  I-l=CONST  LINE 

SL0I11=(CQD1 (I-l, J)-CQD1 (I-l, J-1) ) /DCTJIJ 
CQI11=CQD1 (I-l, J-1) +SL0I11*DCTJ1X 
SLOI12=(CQD2 (I-l, J) -CQD2 (I-l, J-1) ) /DCTJIJ 
CQI12=CQD2 (I-l, J-1)+SL0I12*DCTJ1X 
C  INTERPOL  FOR  CQ1,CQ2  W.R.T.  MU,  CT=CONST 

DMUI 1 I=DMU ( I ) -DMU ( I-l ) 

DMUI 1X=XMU-DMU ( I- 1 ) 

SLCQ1= (CQIl-CQIll) /DMUIII 
CQ1=CQI11+SLCQ1*DMUI1X 
C  INTERPOL  FOR  CQ2  W.R.T.  MU,  CT=CONST 
SLCQ2= (CQI2-CQI12) /DMUIII 
CQ2=CQI 12+SLCQ2  *DMUI IX 
C  CQE5=CQ1+ (CQ2-CQ1) * (SCN-245. 6) /12 . 4 
VTIP1=713.65 
VTIP2=749.61 

CQE5=CQ1+ (CQ2-CQ1) * (VTIP-VTIPl) / (VTIP2-VTIP1) 

CQ=CQE5*1.0E-5 

RETURN 

END 

SUBROUTINE  PAVAIL (NFV, ZNOENG, UTL, XLTL) 
COMMON/PAVA/PAC1E,PAM1E,VF(20) ,PAC1 (20) ,PAM1 (20) ,PAC2 (20) , PAM2 (20) 
DO  100  1=1, NFV 
PACl (I)=PAC1E 
PAMl (I)=PAM1E 

IF(VF(I)  .GE.  50.0)PACl(I)=PACl(I)+0.42*(VF(I)-50.0) 

IF(VF(I)  .GE.  50.0)PAM1  (I)=PAM1  (I)+0.62*  (VF(I)-50.0) 

PAC2 ( I ) =ZNOENG*PACl ( I ) 

PAM2 (I) =ZNOENG*PAMl (I) 

PACl (I)=AMIN1 (PACl (I) , XLTL) 

PAMl (I)=AMIN1 (PAMl (I) ,XLTL) 

PAC2 (I)=AMIN1 (PAC2 (I) ,UTL) 

PAM2 ( I ) =AMIN1 (PAM2 ( I ) , UTL) 

10-  CONTINUE 
RETURN 
END 
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o  o  o 


C 


SUBROUTINE  CTINT 

THIS  SUBRROUTINE  CALCULATES  LOAD  FACTOR, BANK  ANGLE, TURN  RADIUS, AND  RA 
OF  TURN, AND  LATERAL  ACCELERATION  FOR  THE  HORIZENTAL  FLIGHT  MODE  ONLY. 
THIS  SUBROUTINE  INTERPOLATES  ON  MU  AND  CQ  TO  FIND  THE  CT  CORRES.  TO  PA 
COMMON  DMU(15) ,CTD(12) ,CQD1 (14,7) , CQD2 (14,7) ,NMUD,NCTD 
C0MM0N/PAVA/PAC1E,PAM1E, VF (20) ,PAC1 (20) ,PAM1 (20) , PAC2 (20) , PAM2 (20) 
COMMON  /BCT/XMU,PSFAC, SIGMA, VTIP, CT, CLSTL,VFTS, RHOl 
COMMON  /BCT1/CTCP,CTMP, IX, PRHF (20) , PCLF (20) , PMLF (20) ,EFF 
120  FORMAT (1H0,2F1 0.2, 4F12 . 2, 2X, E12 . 4, 3F11.2) 

TSFAC=PSFAC*XMU*550 . 0/VFTS 
WTT=CT*TSFAC 
VFKN=VFTS/1.689 
RCS=0 . 0 

CPCE5=PAC2 (IX) *1 . 0E+5/PSFAC 
CPME5=PAM2 (IX) *1 . 0E+5/PSFAC 
XCP=CPCE5*EFF 

FOR  K=1 , PAV=PAC2 ( IX) , FOR  K=2 , PAV=PAM2 ( IX) . 

DO  90  K=l,2 

IF(K  .EQ.  2)XCP=CPME5*EFF 
10  DO  20  1=1,14 

IF(XMU  .GT.  DMU(I))GOTO  20 
GOTO  30 
20  CONTINUE 
1=1-1 

30  DO  50  J=l,7 

IF(XCP  .GT.  CQDl (I, J) )GOTO  50 
GOTO  60 
50  CONTINUE 
J=J-1 

60  IF(J  .EQ.  1)J=J+1 

SLOPIl= (CTD ( J) -CTD ( J-1) ) / (CQDl (I, J) -CQDl (I, J-1) ) 

CTI1=CTD (J-1) +SLOPI1* (XCP-CQDl (I,  J-1)  ) 

IF(I  .EQ.  1) 1=1+1 
DO  65  JJ=1,7 

IF(XCP  .GT.  CQDl (I-1,JJ) ) GOTO  65 
GOTO  68 
65  CONTINUE 
JJ=JJ-1 

68  IF(JJ  .EQ.  1)JJ=JJ+1 

SLOPIll=(CTD(JJ) -CTD (JJ-1) )/ (CQDl (I-l, JJ) -CQDl (I-l, JJ-1) ) 

CTI11=CTD ( JJ-1) +SLOPI11* (XCP-CQDl (I-l, JJ-1) ) 

CT1=CTI11+(CTI1-CTI11) * (XMU-DMU(I-l) ) / (DMU (I) -DMU (I-l) ) 

DO  70  J=l,7 

IF(XCP  .GT.  CQD2 (I, J) )GOTO  70 
GOTO  75 
70  CONTINUE 
J=J-1 

75  IF(J  .EQ.  1)J=J+1 

SLOP 12= (CTD ( J) -CTD ( J-1) ) / (CQD2 (I, J) -CQD2 (I, J-1) ) 

CTI2=CTD (J-1) +SLOPI2* (XCP-CQD2 (I, J-1) ) 

DO  76  JJ=1,7 

IF(XCP  .GT.  CQD2 (I-l, JJ) ) GOTO  76 
GOTO  77 

76  CONTINUE 
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JJ=JJ-1 

77  IF(JJ  .EQ.  1)JJ=JJ+1 

SL0PI12=(CTD(JJ)-CTD(JJ-1) ) / (CQD2 (I-l, JJ) -CQD2 (I-l, JJ-1) ) 
CTI12=CTD(JJ-l)+SLOPI12* (XCP-CQD2 (I-l,  JJ-1) ) 

CT2=CTI12+ (CTI2-CTI12) * (XMU-DMU (I-l) ) / (DMU (I) -DMU (I-l) ) 

C  CTE4=CT1+ (SCN-245 . 6) * (CT2-CT1) /12 . 4 
VTIP1=713.65 
VTIP2=749.61 

CTE4=CT1+ (CT2-CT1) * (VTIP-VTIPl) / (VTIP2-VTIP1) 

CTX=CTE4*  l.OE-4 
CTXOr  :</ SIGMA 

IF(C  j  .GT.  CLSTL)CTX=CLSTL* SIGMA 
IF(K  .C.Q.  2)  GOTO  80 

C  IN  THIS  PART  WE  CALCULATETHE  LOAD  FACTOR, BANK  ANGLE, TURN  RADIUS, 

C  RATE  OF  TURN  AND  LATERAL  ACCELERATION. 

PCLF (IX)=CTX/CT 

IF(PCLF(IX)  .LT.  1.0)GOTO  81 

TANPHI=SQRT(PCLF(IX)  *PCLF(IX)  -1.0) 

ACCEL=32 . 17*TANPHI 
GOTO  83 

81  TANPHI=-1.0E-5 
ACCEL=32 . 17*TANPHI 

83  PHIR=ATAN(TANPHI) 

PHID*PHIR*57 .29578 
IF(TANPHI  .EQ.  0.0) GOTO  78 
TURRFT=VFKN*VFKN/ (11 . 26*TANPHI) 

GOTO  79 

78  TURRFT=1000000. 

79  ROTDS=1091.*TANPHI/VFKN 

WRITE (*, 120) VFKN,RCS,PRHF (IX) ,PAC2 (IX) ,PCLF(IX) , PHID, TURRET, ROTDS, 
*ACCEL 
GOTO  90 

80  PMLF(IX)=CTX/CT 
IF(PMLF(IX)  .LT.  1.0)  GOTO  86 
TANPHI=SQRT (PMLF (IX) *PMLF (IX) -1.0) 

ACCEL=32 . 17*TANPHI 

GOTO  87 

86  TANPHI=-1 . OE-5 
ACCEL=32 . 17*TANPHI 

87  PHIR=ATAN(TANPHI) 

PHID=PHIR*57 .29578 
IF(TANPHI  .EQ.  0.0) GOTO  85 
TURRFT=VFKN*VFKN/ (11 . 26*TANPHI) 

GOTO  82 

85  TURRFT=1 000000. 

82  ROTDS=1091.0*TANPHI/VFKN 

WRITE(*,120)VFKN,RCS,PRHF(IX) ,PAM2 (IX) , PMLF (IX) , PHID, TURRET, ROTDS, 
*ACCEL 
GOTO  90 
90  CONTINUE 
RETURN 
END 
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SUBROUTINE  VINDRA(NFV,VF,PRHF) 

DIMENSION  XD(5)  ,X(5,3)  ,  Y(5)  ,  AF  (5)  ,  A(4, 4)  ,  C  (3)  ,SIG(3)  ,T(3)  ,R(5)  , 
*XC(2)  ,  YC(2)  ,VF(1)  ,PRHF(1) 

C  THIS  SUBROUTINE  CALCULATES  THE  INDURANCE  VELOCITY  AND  THE  VELOCITY 

C  CORRESPONDING  TO  MAXIMUM  RANGE. THE  COMPUTATION  IS  PERFORMED  WITH  A 

C  LEAST  SQUARE  FIT  OF  FIVE  DATA  POINTS  OF  THE  POWER  VS  VELOCITY  CURVE 

IC=0 
N=3 
M=5 
NRX=5 
NRA=4 
IND=NFV-1 
DO  50  J=1,IND 

IF(PRHF(J-*-l)  .LT.  PRHF(J))GOTO  50 
GOTO  60 
50  CONTINUE 
60  XD(l)=VF(J-2) 

XD(2)=VF(J-1) 

XD(3)=VF(J) 

XD(4)=VF(J+1) 

XD(5)=VF(J+2) 

Y(l)=PRHF(J-2) 

Y(2)=PRHF(J-1) 

Y(3)=PRHF(J) 

Y(4)=PRHF(J+1) 

Y(5)=PRHF(J+2) 

DO  100  1=1, M 

X(I,1)=1 

X(I,2)=XD(I) 

X{I,3)=XD(I)*XD(I) 

100  CONTINUE 

CALL  GENLSQ (X, NRX, Y, M, A, NRA, N, C, R, AF, ERMS, SIG, T, DET, IC) 

WRITE(*,2) 

WRITE(*,1)C 

WRITE(*,2) 

WRITE(*,3)  (XD(I)  ,Y(I)  ,AF(I),R(I),I=1,M) 

WRITE(*,2) 

XC(1)=-C(2)/(2.0*C(3)) 

XM=C(2)+2,0*SQRT(C(1)*C(3)  ) 

XC(2)  =  (XM-C(2) )  /  (2.0*C(3) ) 

DO  200  K=l,2 

200  YC  (K)  =C  (1)  +C  (2)  *XC  (K)  +C  (3)  *XC  (K)  *XC  (K) 

1  FORMATdHO,  'C(l)='  ,F10.5,5X,  'C(2)=',F10.5,5X,  'C(3)='  ,F10.5) 

2  FORMAT (//) 

3  FORMAT (15X,4F1 0.5) 

4  FORMAT (IHO, ' VMAXIND=' , FIO . 5, 5X, ' PRMIN=' , FIO . 5, 5X, ' VMAXRA=' , FIO . 5, 
*5X, ' PRMAXRA=S FIO . 5) 

WRITE(*,4)  (XC(I)  ,YC(I)  ,1=1,2) 

RETURN 

END 

SUBROUTINE  LACC (CPCE5, CPME5, VTIP, CTCP, CTMP) 

COMMON  DMU(15) ,CTD(12) ,CQD1 (14,7) ,CQD2 (14,7) ,NMUD,NCTD 
C  THIS  SUBROUTINE  CALCULATES  THE  MAX  AMOUNT  OF  THRUST  THAT  CAN 
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BE  GENERATED  IN  HOVER  FOR  THE  POWERS  PAC,AND  PAM. 

CQX=CPCE5 
0  90  K=l,2 

IF(K  .EQ.  2)CQX=CPME5 
DO  50  J=l,7 

IF(CQX  .GT.  CQDl (1, J) )GOTO  50 
GOTO  60 
50  CONTINUE 
60  IF(J  .EQ.  1)J=J+1 

CT1=CTD (J-1) + (CTD (J) -CTD (J-1) ) * (CQX-CQDl (1, J-1) ) / (CQDl (1, J) -CQDl (1 
J-1)  ) 

DO  70  JJ=1,7 

IF(CQX  .GT.  CQD2 (1, JJ) )GOTO  70 
GOTO  80 
70  CONTINUE 

80  IF(JJ  .EQ.  1)JJ=JJ+1 

CT2=CTD ( JJ- 1 )  +  (CTD ( JJ) -CTD ( JJ-1 ) ) * (CQX-CQD2 ( 1 ,  JJ-1 ) ) / (CQD2 ( 1 ,  JJ) -C 
*QD2 (1, JJ-1) ) 

CTXE4=CT1+(CT2-CT1) * (SCN-245.6) /12.4 
VTIP1=713.65 
VTIP2=749.61 

CTXE4=CT1+ (CT2-CT1) * (VTIP-VTIPl) / (VTIP2-VTIP1 ) 

CTX=CTXE4*1.0E-4 
IF(K  .EQ.  1)CTCP=CTX 
IF(K  .EQ.  2)CTMP=CTX 
90  CONTINUE 
RETURN 
END 


SUBROUTINE  GE (CT, CQI , FL, CPCE5, CPME5 , CTCI , CTMI ) 

COMMON/IGE/CTIN(10) ,CQIN(10) ,XINMU(10) 

IF(FL  .EQ.  1.0) GOTO  40 

THIS  SUBROUTINE  CALCULATES  IN  THIS  PART  THE  CQ  REQUIRED  FOR  HOVER  IN  I 

CTE4=CT*1.0E+4 

DO  20  1=1,8 

IF(CTE4  .GT.  CTIN(I))GOTO  20 
GOTO  30 
20  CONTINUE 
1=1-1 

30  IF(I  .EQ.  1) 1=1+1 

CQE5=CQIN (I-l) + (CQIN (I) -CQIN (I-l) ) * (CTE4-CTIN (I-l) ) / (CTIN (I) -CTIN ( 
*1-1)) 

CQI=CQE5*1.0E-5 
GOTO  100 

THIS  SUBROUTINE  CALCULATES  IN  THIS  PART  THE  THRUST  IT  CAN  GENERATE  IN 
THE  I.G.E.  CORRESPONDING  TO  PAC  AND  PAM. 

40  CQX=CPCE5 
DO  90  K=l,2 
IF(K  .EQ.  2)CQX=CPME5 
DO  60  J=l,8 

IF(CQX  .GT.  CQIN (J)) GOTO  60 
GOTO  70 
60  CONTINUE 


70  IF(J  .EQ.  1)J=J+1 

CTX4*CTIN(J-1)+ (CTIN(J)-CTIN(J-l) ) * (CQX-CQIN ( J-1 ) ) / (CQIN(J) -CQIN(J 
*-l)) 

CTX=CTX4*1.0E-4 
IF(K  .EQ.  1)CTCI=CTX 
IF(K  .EQ.  2)CTMI=CTX 
90  CONTINUE 
100  RETURN 
END 

SUBROUTINE  GENLSQ (X, NRX, F, M, A, NRA, N, C, R, AF, ERMS, SIG, T, DET, IC) 
DIMENSION  X(NRX,N)  ,F(NRX)  ,  A  (NRA,  N)  ,  C  (N)  ,R(M)  ,AF(M)  ,SIG(N)  ,T(N)  , 
ISA(IOO) 

NN=N 
MM=M 
NP1=NN+1 
DO  9  IJ=1,NN 
DO  9  IK=1,NP1 
A(IJ,IK)=0.0 
RSQAR=0 . 

10  DO  12  K=1,MM 
J=1 

DO  11  1=1, NN 
SA(J)  =X(K,  I) 

11  J=J+1 
SA(J)=F(K) 

W=l. 

12  CALL  FNEQS(A,NN,SA,NRA,W) 

CALL  MATINV(A,NN,C,NRA,2,DET) 

DO  125  1=1, NN 

125  C(I)=A(I,NN+1) 

IFdC.EQ.DGO  TO  18 
DO  14  K=1,MM 
APP=0. 

DO  13  1=1, NN 

13  APP=APP+X(K, I) *C(I) 

TE=F (K) “APP 
RSQAR=RSQAR+TE * TE 
IF(IC.EQ.2)GO  TO  14 
R(K)=TE 

IF(IC.EQ.3)GO  TO  14 
AF  (K)  =APP 

14  CONTINUE 
IF(M.LE.N)  GO  TO  19 
FM»MM 
FN=NN 

ERMS=SQRT  (RSQAR/  (FM-FN)  ) 

DO  17  1=1, NN 
IF(A(I,I)  .LE.0.)GO  TO  16 
SIG(I)=ERMS*SQRT(A(I,I)  ) 

T(I)=C(I)  /SIG(I) 

GO  TO  17 
16  SIG(I)=A(I,I) 

T(I)=0. 
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17  CONTINUE 

18  RETURN 

19  PRINT  20,MM,NN 
ERMS=0 . 

DO  30  1=1, NN 
SIG(I)=0. 

30  T(I)=0. 

GOTO  18 

20  FORMAT (23H  ERROR-GENLSQ  M.LE.N  M=, 15, 3H  N=, 15) 
END 

SUBROUTINE  MATINV (A, N, C, NMAX, K, DET) 

DIMENSION  A (NMAX,N) ,  C (N) 

NN=N 
KK=K 
ZERO=0 . 

ASSIGN  30  TO  N6 
IF(l-KK)  24,23,23 

23  N3=NN 

IF(KK)  25,28,25 

28  ASSIGN  40  TO  N6 
GO  TO  26 

25  ASSIGN  8  TO  N5 
ASSIGN  16  TO  N7 
ASSIGN  10  TO  N8 
GO  TO  27 

24  N3=KK+NN-1 
ASSIGN  2  TO  N8 

26  ASSIGN  9  TO  N5 
ASSIGN  17  TO  N7 

27  T3=1.0 
T2=1.0 

DO  21  1=1, NN 

IF(ABS(A(I,I)  )-ZERO)  29,29,2 

29  GO  TO  N6,  (40,30) 

40  PRINT  4 

DET  =0.0 
RETURN 

4  FORMAT (16H  SINGULAR  MATRIX) 

2  T1=1.0/A(I,I) 

T3=T3*A(I, I) 

A(I,I)=1.0 
DO  6  J=1,N3 

6  A(I,  J)=A(I,  J)  *T1 

7  GO  TO  N5,  (8,9) 

8  C(I)=C(I)*T1 

9  DO  17  J=1,NN 
IF(I-J)  11,17,11 

11  T1*A(J,  I) 

A(J,I)=0.0 
DO  14  L=1,N3 

14  A(J,L)=A(J,L)-T1*A(I,L) 

15  GO  TO  N7,  (16,17) 

16  C(J)=C(J)-T1*C(I) 
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17  CONTINUE 

21  CONTINUE 
IF(T3)  22,50,22 

50  T3=1.0E-150 

22  DET=T3*T2 
RETURN 

30  IF(I-NN)  31,40,31 

31  J=I+1 

DO  39  L=J,NN 

IF(ABS(A(L,I)  )-ZERO)  39,  39,42 
39  CONTINUE 
GO  TO  40 
42  DO  41  J=1,N3 
T1=A(I,  J) 

A(I,  J)=A(L,  J) 

41  A(L,J)=T1 
T2=-T2 

GO  TO  N8,  (10,2) 

10  T1=C(I) 

C(I)=C(L) 

C(L)=T1 
GO  TO  2 
END 

SUBROUTINE  FNEQS (A,N,C,NMAX,W) 
DIMENSION  A(NMAX,N)  ,C(N) 

JA=N+1 

L=1 

DO  10  1=1, N 

M=L 

K=I 

T=  C(I)*W 
IF  (T)5,10,5 
5  DO  14  J=1,JA 
IF  (J-I)9,8,8 

8  A(K,1)=A(K,1)+T*C(J) 

14  K=  K+NMAX 

10  L=  L+NMAX 
RETURN 

9  A(K  ,1)=A(M,1) 

M=M+1 

GOT014 

END 
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APPENDIX  B; 

USTINO  OF  TORQUE  VS.  THRUST  COEFHOENTS  OF  TEST  DATA 
FOR  THE  UH-60  AND  THE  AH64  HEUCOPTERS 


INTENTIONALLY  LEFT  BLANK. 
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UH60Data 

OUT  OF  GROUND  EFFECT  TORQUE  COEFFICIENTS, 
CQ  *  10*  vs.  CT  *  10^  and  ji 


1  NUMUD=  14 

NCTD  =  7  1 

Ii/cr 

66.0 

72.0 

78.0 

mm 

96.0 

102.0  1 

0.00 

63.2 

70.9 

79.0 

87.3 

95.9 

104.9 

114.1 

0.12 

34.5 

38.5 

43J 

48.4 

54.0 

60.2 

68.2 

0.14 

34.0 

37.2 

41.5 

46.0 

50.7 

56.5 

63.5 

0.16 

33.0 

36.2 

39.7 

43.8 

48.2 

53J 

60.0 

0.18 

32.6 

35.5 

38.5 

42.5 

46.7 

51.5 

57.8 

0.20 

32.7 

35.2 

38.4 

42.4 

46.5 

51.4 

57.1 

0.22 

33.5 

36.0 

39.2 

43.0 

412 

52.0 

58.0 

0.24 

35.2 

37.8 

41.0 

45.0 

492 

54.2 

60.5 

0.26 

38.3 

40.7 

44.0 

48.0 

52.6 

58.0 

65J 

0J8 

41.7 

44J 

48.5 

52.5 

57.7 

64.2 

73.0 

0.30 

47.0 

49.6 

54.0 

59.0 

65.0 

73.0 

86.0 

0.32 

54.0 

56.1 

61.0 

66.3 

142 

85.0 

102.5 

0.34 

62.0 

64.2 

69J 

76.1 

85.3 

99.0 

130.0 

0.36 

71J 

74.0 

80.0 

87.5 

100.1 

118.5 

155.0 
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UH60  Data  (continued) 

OUT  OF  GROUND  EFFECT.  SECOND  SET  OF  TORQUE  COEPHCIENTS, 

CQ  *  10^  vs.  CT  •  10* 


H/CT 

66.0 

72.0 

78.0 

84.0 

90.0 

96.0 

102.0 

0.00 

64.6 

72.6 

80.9 

89.5 

98.5 

107.7 

117.3 

0.12 

35.0 

38.5 

43.5 

49.0 

54.5 

61.5 

70.5 

0.14 

34.2 

37.5 

41.7 

45.7 

50.5 

56.0 

63.0 

0.16 

33.0 

36.0 

39.5 

43.5 

48.0 

53.5 

59.0 

0.18 

32.5 

35.5 

38.5 

42.5 

46.5 

51.0 

56.0 

0.20 

32.6 

35.6 

39.0 

41.8 

45.5 

50.0 

54.7 

0.22 

33.5 

35.8 

39.0 

42.2 

36.0 

50.2 

55.0 

0.24 

35.2 

37.7 

40.5 

43.8 

47.6 

52.0 

57.5 

0.26 

37.7 

40.2 

43.0 

46.5 

51.0 

56.0 

62.0 

0.28 

41.1 

43.5 

46.5 

50.5 

55.5 

62.0 

69.5 

0.30 

46.0 

48.0 

51.0 

55.5 

62.0 

70.0 

80.0 

0.32 

52.0 

54.0 

515 

63.0 

70.5 

80.5 

92.5 

0.34 

59.7 

61.5 

65.5 

72.0 

81.0 

95.0 

112.5 

0.36 

68.8 

71.0 

75.5 

82.5 

94.0 

115.0 

130.0 

UH60  Data  (continued) 

INGROUND  EFFECTS  TORQUE  COEFFICIENTS, 
CT  *  10*  vs.  CQ  *  10^ 


CT*  10* 


60.0 


66.0 


72. 


78. 


90.0 


96.0 


102.0 


CQ*  10* 


45.0 

0.0 

50.2 

0.0 

56.0 

0.0 

62J 

0.0 

165 

0.0 

84.0 

0.0 

92.5 

0.0 
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AH64Data 

OUT  OF  GROUND  EFFECT  TORQUE  COEFHCIENTS. 
CQ  *  10^  vs.  CT  *  10*  AND 


1  NMUD=  14 

NCTD  =  7  1 

1  ^ 

50.0 

59.0 

68.0 

77.0 

86.0 

95.0 

104.0 

0.00 

46.0 

58.0 

71.1 

84.6 

98.8 

113.8 

129.6 

0.12 

32.3 

36.2 

41.4 

48.0 

56.0 

64.5 

74.9 

0.14 

30.0 

33.9 

38.8 

45.1 

52.6 

60.8 

70.2 

0.16 

28.4 

32.3 

37.1 

433 

50.8 

58.7 

67.8 

0.18 

27.6 

31.5 

36.3 

42.5 

49.6 

57.5 

66.7 

0.20 

30.0 

33.8 

38.1 

43.6 

50.1 

57.9 

67.9 

0.22 

32.4 

35.9 

40.4 

46.1 

52.6 

60.4 

70.2 

0.24 

34.9 

38.7 

43.1 

48.6 

553 

63.7 

74.5 

0.26 

39.6 

43.0 

47.2 

52.9 

59.7 

68.3 

79.8 

0.28 

44.8 

48.2 

52.9 

58.7 

65.7 

75.0 

88.3 

0.30 

51.8 

55.8 

60.5 

66.2 

73.6 

83.3 

99.6 

0.32 

60.6 

64.9 

70.2 

77.0 

85.5 

96.4 

116.4 

0.34 

71.3 

76.1 

82.0 

89.4 

99.0 

112.4 

137.0 

0.36 

81.6 

89.0 

95.9 

104.9 

116.9 

133.7 

163.3 
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AH64  Data  (continued) 

OUT  OF  (MOUND  EFFECT,  SECOND  SET  OF  TORQUE  COEFHOENTS, 

CQ  •  IC?  vs.  CT  *  10* 


p/CT 

50.0 

59.0 

68.0 

77.0 

86.0 

95.0 

104.0 

0.00 

46.9 

58.5 

71.1 

84.6 

98.8 

113.8 

129.6 

0.12 

32.3 

36.2 

41.4 

48.0 

56.0 

64.5 

74.9 

0.14 

30.0 

33.9 

38.8 

45.1 

52.6 

60.8 

70.2 

0.16 

28.4 

32.3 

37.1 

43.5 

50.8 

58.7 

67.8 

0.18 

27.6 

31.5 

36.3 

42.5 

49.6 

57.5 

66.7 

0.20 

30.0 

33.8 

38.1 

43.6 

50.1 

57.9 

67.9 

0.22 

32.4 

35.9 

40.4 

46.1 

52.6 

60.4 

70.2 

0.24 

34.9 

38.7 

43.1 

48.6 

55.5 

63.7 

74.5 

0.20 

39.6 

43.0 

47.2 

52.9 

59.7 

68.3 

79.8 

0.28 

44.8 

48.2 

52.9 

58.7 

65.7 

75.0 

88.3 

0.30 

51.8 

55.8 

60.5 

66.2 

73.6 

83.3 

99.6 

0.32 

60.6 

64.9 

70.2 

77.0 

85.5 

96.4 

116.4 

0.34 

71.3 

76.1 

82.0 

89.4 

99.0 

112.4 

137.0 

0.36 

81.6 

89.0 

95.9 

104.9 

116.9 

133.7 

163.3 

AH64  Data  (cratinued) 

INGROUND  EFFECTS  TORQUE  COEFFDBIENTS, 
CT  ♦  10*  vs.  CQ  *  10* 


Cr*  10* 

CQ*  10* 

'  ! 

41.0 

30.7 

0.0 

1 

39.4 

0.0 

59.0 

49.0 

0.0 

68.0 

59.3 

0.0 

77.0 

70.3 

0.0 

86.0 

82.0 

0.0 

95.0 

94.3 

0.0 

104.0 

1072 

0.0 

46 


UST  OF  SYMBOLS 


VARIABLE 

HELPE  PROGRAM 

SYMBOL 

NAME 

h 

Operating  altitude  in  feet 

ALTl 

IRP 

Intennediate  rated  power  in  horsepower. 

TIRP.PAM 

MCP 

Maximum  continuous  power  in  horsepower. 

TMCP,  PAC 

Vk 

Fbrward  speed  (knots). 

V 

Fbrward  speed  (ft/isec). 

VFTS 

Vt 

Rotor  tip  speed  in  ft/sec. 

VTIP 

1* 

Ratio  of  forward  qreed  to  rotor  tip  speed  (ND). 

W.  DMU,  XMU 

W 

Operating  weight  0n  pounds). 

WT 

R 

Main  rotor  radius  (feet). 

ROTR 

N 

Number  of  rotor  revolutions  per  second. 

N 

Po 

Power  loss  due  to  blade’s  profile  drag. 

A 

Rotor  disk  area  (f^). 

AD 

Af 

Frcmtal  drag  area  of  a  helicopter  (ft^). 

FA 

At 

Tq)  flat  plate  area  of  a  helicopter  (fr). 

TFPA 

Cd 

Coefficiem  of  drag. 

CD,  CDX 

f 

Equivalent  fiontal  drag  area  (C^  •  Af). 

F 

P 

Air  den^. 

RHOI 

CQ 

Coeffident  of  torque  (ND). 

DCQ 

Incremental  torque  coefficient  due  to  drag. 

DCQ 

CP 

Coefficient  of  power  (ND). 

CT 

Coefficient  of  dmist  (ND). 

CT 

Psc 

Power  scale  factor  fibs  ft^). 

PSFAC 

Tsc 

Thrust  scale  factor  flbs). 

TSFAC 

Ph 

Power  required  for  horizmital  flight  (HP). 

PRHF 

Pc 

Pow^  required  for  vertical  dimb  (HP). 

PRCL 

Pt 

Total  power  required  (HP). 

PRTOT 

If 

Load  factor  0^). 

XKF1,XLF2 

♦ 

Bank  angle  in  a  turn. 

PHIr.PHID 

T 

Rotor  ditust  flbs). 

Rt 

Turn  racUus  (ft). 

TURRFT 

6 

Rate  of  turn  in  deg/sec. 

ROTDS 

Pc 

Centripetal  force  (ft^sec?). 

Number  of  missimi  wd^its  to  compute. 

NWT 

Number  of  levd  flight  qreed  considered. 

NFV 

Number  of  vertical  clunb  rates. 

NCR 

Levd  speed  array  name  (knots). 

VF 

PACl 

Qmtinuous  engine  power  avaMile  (HP). 

PACl 

PAMl 

Maximum  engine  power  available  (Iff*)- 

PAMl 

e 

Rotor  system  effidency  foctor  (ND). 

EFF 

Rate  of  dimb  (sec*^). 

RCS 
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